Atomistic approach to alloy scattering in Sii-xGe. 
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SiGo alloy scattering is of significant importance with the introduction of strained layers and SiGe 
channels into CMOS technology. However, alloy scattering has till now been treated in an empirical 
fashion with a fitting parameter. We present a theoretical model within the atomistic tight-binding 
representation for treating alloy scattering in SiGe. This approach puts the scattering model on a 
solid atomistic footing with physical insights. The approach is shown to inherently capture the bulk 
alloy scattering potential parameters for both n-type and p-type carriers and matches experimental 
mobility data. 



The SiGe alloy has found widespread application 
in complementary metal-oxide semiconductor (CMOS) 
technology [1]. Continued scaling of CMOS devices has 
led to the exploration of new materials to augment the 
silicon channel. Considerable effort has been spent on 
the realization of ultrathin strained SiGe channel devices 
[2, 3] with alloy scattering being one of the major scat- 
tering mechanisms limiting the mobility [4, 5]. 

Alloy scattering has previously been modeled as scat- 
tering due to random potential wells of radius R and 
depth AC/ [6]. The uncertainity in the definition and 
large discrepancy in the reported values of AU has widely 
been discussed in literature [4, 5]. Due to the lack of a 
consistent physical basis, AU has been relegated to a 
mere fitting parameter. The work of Fischetti and Laux 
[4] is noteworthy in providing a set of experimentally fit- 
ted AU values. There have been some recent theoreti- 
cal effort in estimating scattering potentials in bulk SiGe 
using first-principles calculations which render transport 
calculations computationally burdensome [7, 8]. 

Yet further critical questions remain to be tackled. How 
do alloy scattering rates vary with strain? How does 
quantum confinement in a system with countable num- 
ber of atoms affect mobility due to alloy scattering? Can 
the fitting parameter AU be filled with a physical mean- 
ing and insight? To address these questions we consider 
a fresh look at alloy scattering from an atomistic per- 
spective. 

We study the treatment of alloy scattering in SiGe us- 
ing a 20 orbital basis set sp^d^s* semi-empirical tight- 
binding (TB) approach [9]. Previously we developed 
a virtual crystal approximation (VGA) bandstructure 
model to study transport in SiGe [10]. The model was 
benchmarked against experimental data for bulk band 
edges and band gaps. Here wc analyse an atomistic al- 
loy disorder representation as a perturbation over the 
homogenous VGA potential. A VGA TB Hamiltonian 
for a Sii-xGcx composition considers "averaged" atomic 
species. Each bond length in VGA is assumed to be equal 
to the averaged Si and Ge bond length for the corre- 
sponding composition x (also known as 'Vegard limit'). 
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FIG. 1. (a) One of the 32 {Nc = 2^) possible Sii-^Gex unit 
cells in an atomistic representation occuring with probability 

p{x) = (1 — x)^x^ . (b) Unit cell representation in VGA with 
all bonds and atoms as averaged quantities. 

As opposed to VGA an atomistic TB Hamiltonian repre- 
sents each 'atom type' as either Si or Gc and each 'bond 
type' as either Si-Si, Si-Ge or Ge-Ge. In a realistic bulk 
relaxed SiGe alloy each bond length tries to assume its 
original bulk bondlcngth (also known as 'Pauling limit') 
[11]. In our calculations, wc take each bond to be equal 
to its bulk bond length while Si-Ge bonds are taken to 
be half the pure Si and Ge bond lengths. Electronic cou- 
pling parameters for the Si-Ge bonds are set to be equal 
to that of VGA S'zo.sGeo.s. 

The bulk atomistic Hamiltonian can now be repre- 
sented as a sum of the VGA Hamiltonian and a disorder 
term. Within first order perturbation theory the disorder 
term, H"^'"'^ exactly contributes to the alloy scattering. 
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The elastic alloy scattering rate t^{E) according to 
Fermi's Golden rule for carriers from a band edge a at 
wave vector fc„ to a band edge (5 at wave vector k'^ can 
now be written as, 

SE(k^)M^'p)^\ka.lW^\ (2) 

where N is the total number of atoms in the bulk Hamil- 
tonian. The TB wavefunction at band edge a is given 
by 
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TABLE I. Extracted alloy scattering potentials AUtb{x) at the CB and VB edges and compared with experimental values. 
The last column hsts the relative change in AUtb with ofT-diagonal block in Eq. 5 set to zero ( \AUTB-AUTB(offdiagonai=o)\ j^ 
where D,0, and B stand for diagonal, off-diagonal and both. AUtb is the best constant fit of AUtb{x). 



Band edge (Sii-xGex) 


Alloy Scattering potential (eV) 

AUtb{x) AUtb AUg^pt/theory* 


Main Contributing Block 
(relative change) 


CB A (xjO.84) 


Omx^- 0.181a;+0.963 


0.90 


0.7 


B(35%) 


CB L (a; > 0.84) 


-0.05a;^-0.08a;+1.78 


1.67 


1.0,0.6 


0(75%) 


CB A - L{x w 0.84) 


0.914 


0.914 






VB r point 


0.777 - 0.06a; 


0.75 


0.9,0.81* 


D(12%) 



where (j){R — r„) represents the orthogonalized Lowdin 
orbital at site r„ (a, /3 G A or L for the conduction band 
(CB) and T for the valence band (VB)). We assume the 
case of a perfect random alloy which leads to the condi- 
tion Eq. 4 [8]. 
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Using Eq. 3 and 4, Eq. 2 can be written in a condensed 
form as the sum of scattering due to Nc (=32) different 
bulk atomistic unit cell types (Figure 1 (a)). 
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Here x is the ratio of Ge content, fio (= a^/S) is the 
unit cell volume where ao is the Sii-xGcx lattice con- 
stant in Vcgard limit. Ny is the number of degenerate 
valleys which is 6 for A, 4 for L and 1 each for heavy and 

light hole band edge. hI (= - HX^^) is the differ- 
ence between the j*'* atomistic unitcell Hamiltonian, Hq 
that occurs with probability pj (x) and the VGA unitcell 
Hamiltonian, H^'-''^. ipo{ka) is the unit cell wavefunc- 
tion for band edge a and Pf^{E) is the density of states 
per spin of the band edge j3 where carriers scatter into 
(q:=/3 for intravalley scattering while a ^ /3 for interval- 
ley scattering). Eq. 5 can be written in an expression 
(Eq. 6) that is commonly used in the literature for alloy 
scattering rates [4]. 
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The effective alloy scattering potential AUtb{x) em- 
bedded in the atomistic Hamiltonian can now be calcu- 
lated using Eq. 5 and 6 as shown in Eq. 7. 
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The extracted potentials AUtb{x) a-re compared with 
experimentally fitted values in Table I. The Scattering 



potential for VB is averaged over all bands at F. A single 
A — L scattering potential is quoted which has been av- 
eraged between A and L density of states. In the view of 
the empirical nature of the Af7/it,the theoretically com- 
puted AUtb values are found to be reasonably close. 

Scattering potentials for different intervalley and in- 
travalley transitions arc shown in Figure 2 (a) . It is seen 
that / — type scattering in the A valley is non-zero as op- 
posed to what was assumed in previous calculations [4]. 
Also it is seen that inter-valley scattering is nearly as 
strong as intra-valley scattering for L minima. Figure 2 
(b) shows the contributions from different unit cell for A 
intravalley scattering in bulk <S'io.i6Geo.84- The scatter- 
ing potentials are found to depend on both, the atomic 
species and their neighbourhood. This is in contrast to 
the standard model which assumes a hard wall potential 
at each atomic site. Now the understanding of alloy scat- 
tering can be redefined as scattering due to variation in 
'atom types' and 'bond types'. From the last column in 




32 unitcell configurations in Si^ .i^Ge^ ^ 



FIG. 2. (a) Alloy scattering potentials calculated for 
Sii-xGcx {x ^ 0.84) (b) Calculated transition Matrix 
elements \{<l'o{kc,)\Hi - HX'^^\'^o{kc))\ [a € A[ioo]) in 
Sio.i6Geo.84 due to the 32 possible unit cell configurations 
contributing to A intravalley scattering. 
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FIG. 3. Calculated electron and hole mobility at 300K in bulk 
relaxed SiGe. Phonon mobility from [4]. Experimental data 
from ▲ [12] ★ [13] and Q [14]. 



Table I it can be seen that p-type carriers are scattered 
more due to the variation in 'atom type' (Diagonal block) 
while n-type carriers are scattered more due to variation 



in 'bond type' (OfT-diagonal block). 

Our approach to scattering rates is then used to com- 
pute bulk mobility for n-type and p-type carriers in 
Sii-xGcx and compared to experimental data. We use 
the Boltzmann transport equation within the isotropic 
parabolic band approximation similar to [8] . Eq. 8 rep- 
resents the mobility expression for carriers at a conduc- 
tion band edge E" (limits reversed for VB edge). 
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The conductivity {rn^'~'^) and density of states [ra^^^^) 
masses are extracted from the VGA bandstructure [15]. 
The final mobility value due to alloy scattering can be cal- 
culated using y^°-'''-°y = J2a''^afJ'a/n where is the car- 
rier density at the band edge a and n is the total carrier 
density.The bulk SiGe phonon limited mobility ^p'*°"°" 
is taken from [4]. Figures 3 (a) and (b) show the com- 
puted bulk mobility for electrons and holes respectively. 
The final mobility is calculated using Matthiessen's rule, 
1/H = l/^P''°»o" + A close match is found with 

the available experimental mobility data for both elec- 
trons and holes. 

To summarize, we have presented an atomistic ap- 
proach in a tight-binding representation to treat alloy 
scattering in SiGe and verified it for bulk against exper- 
imental data. The theory predicts AU values which are 
close to the experimentally fitted numbers. The bulk mo- 
bility is computed for both n-type and p-type SiGe mate- 
rials using the extracted potentials which closely match 
experimental data. These results indicate the importance 
of atomistic treatment of alloy scattering in SiGe and 
further allow studying its effects in nanostructures in a 
realistic fashion. 
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